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We propose and verify a wave- vector-space version of gen- 
eralized extended self similarity Q and broaden its appli- 
cability to uncover intriguing, universal scaling in the far 
dissipation range by computing high-order (< 20) struc- 
ture functions numerically for: (1) the three-dimensional, in- 
compressible Navier Stokes equation (with and without hy- 
perviscosity); and (2) the GOY shell model for turbulence. 
Also, in case (2), with Taylor-microscale Reynolds numbers 
4x 10'' < Rex < 3x 10", we find that the inertial-range expo- 
nents (^p) of the order - p structure functions do not approach 
their Kolmogorov value p/3 as Re\ increases. 
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The central concern of studies of homogeneous, 
isotropic turbulence is the scaling of order-p velocity 
structure functions, e.g., Sp{r) = (|vi(x + r) — Vj(x)|^'), 
where z(= 1, 2, or 3) is the Cartesian component of the ve- 
locity v(x) at point x, and the angular brackets imply, in 
principle, a spatiotemporal average. Kolmogorov (K41) 
1^ predicted that, at high Reynolds numbers Rex and for 
the inertial range 20?]^ ^ r <^ L (rjd and L are, respec- 
tively, dissipation and forcing scales and A is the Taylor 
microscale), Sp{r) ~ r^p with C,p — p/3. Subsequent ex- 
perimental and theoretical studies [[ip-[ll[| have argued 
for: (1) multiscaling, i.e., Cp — p/3 — S(p, with S^p > 
but Cp a nonlinear, monotonically increasing function of 
p; and (2) extended self similarity (ESS) [^, in which 

is obtained from Sp S^'' , since this extends the ap- 
parent inertial range down to r ~ 5?7rf. A recent gener- 
alization uses Gp{r) = Sp{r)/[S3{r)]P^^ and suggests 
that a log-log plot of Qp versus Qq is a straight line with 
slope pp^q = [Cp—pCa/Sl/iCq — 9(3/3] for the lowest resolv- 
able values of r. This generalized extended self similarity 
(GESS) has been tested ^ to some extent (p, q <G). 

Here we show how GESS is modified at suffi- 
ciently small r by computing wave- vector-space (A:-space) 
analogs of high-order (< 20) structure functions for 

(1) the three-dimensional, incompressible Navier Stokes 
equation {3d NS), with and without hyperviscosity, and 

(2) the GOY shell model for turbulence (where 
we attain both large Re\ and k ^ kd = Vd^)- We fur- 
ther propose a fc-space GESS §, show that it holds for 
L^^ fc ^ 1.5kd, but then crosses over to another form 
in the far dissipation range. To study this we postulate 
/c-space ESS (for real-space structure functions we use 
the symbols S and Q and for their fc— space analogs (not 



Fourier transforms) the symbols S and G): 

Sp = (|v(k)|f ) « AjpiSsf'p, « A: < 1.5kd, 

Sp = {W{k)\P)^ADp{S3r-, 1.5fcd<fc«A, (1) 

where Ajp and Aop are, respectively, nonuniversal am- 
plitudes for inertial and dissipation ranges and A~^ the 
(molecular) length at which hydrodynamics fails (see 
for real-space analogs). Our study shows (Figs|^||) that 
Eq. (1) holds with two different exponents ap and t^p. In 
the GOY model Cp = Cp; but we find explicitly [inset (b), 
Fig.l] that, for the M NS case, Cp = 2(Cp + 3p/2)/ll 
(i.e., Sp{k) - fc-(Cp+3p/2) in i-ijg inertial range @); 
the difference between the two arises because of phase- 
space factors. Both (p and ap (Figj2|) seem universal 
(the same for all GOY and 3d NS runs (Table |) @). 
Cp agrees fairly with the She-Leveque (SL) [Q formula 
(SL ^ 2[1 - (2/3)P/3] for the ranges of p and Rex 

in Fig.||; and ap is close to, but systematically less than, 
p/3. 

The k dependences of the inertial- and dissipation- 
range asymptotic behaviors follow now from the depen- 
dence of 5*3 on k: We find 



^3 ~ Bik-^'-'^/^, L-^ ^k< l.bkd, 
S3 w Bok^ exp{-ck/kd), 1.5kd < A; < A, 



(2) 
(3) 



where Bj and B]j are, respectively, nonuniversal ampli- 
tudes (Eq. (2) holds jl| for 3d NS; for GOY the factor 
9/2 is absent). Thus, in the far dissipation range, all Sp ^ 
k^p exp{—capk/kd) for 1.5kd ^ k A, with 9p — apS, a 
form not easy to verify numerically for large p, given the 
rapid decay at large k, and suggested hitherto jl^ only 
for In Eq. (3), S,c,kd are not universal, but we ex- 
tract the universal part of the crossover via our /c-space 
GESS: Define Gp = Sp/iS3)P/^; log-log plots of Gp versus 
Gq now yield curves (Figs. ^ and ^d) with asymptotes 
which have universal, but different, slopes in inertial and 
dissipation ranges. The inertial-range asymptote has a 
slope p(p, q) (as in real-space GESS § which follows from 
the formulae above); the resulting (^p are in fair agreement 
with the SL formula 0|. The dissipation-range asymp- 
tote has a slope a;(p, g) = [ap—p/3]/[aq~q/3] (see Eq.(l) 
and the definition of Gp). The slopes of these asymp- 
totes are universal, but the point at which the curve 
veers off from the inertial-range asymptote depends on 
the model (GOY, NS, etc.). However, a simple transfor- 
mation yields a universal crossover scaling function (dif- 
ferent for each (p, q) pair because of multiscaling): Define 
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FIG. 1. Log-log plots (base 10) of Sp versus S-i for 3d NS 
(p=17 for runs NSl-4) and GOY (run Gl in inset (a)) models 
showing our fc-space ESS (Eq. 1); full lines are the SL pre- 
diction [^]. Inset (b); (^!p (circles) from run NS4; the line is 
Cp = 2(Cp + 3p/2)/ll, with the (p = Cp^- Note the devia- 
tion of our data points from SL lines at small S3, i.e., in the 
dissipation range. 



FIG. 2. Inertial- and dissipation-range exponents and 
Qfp (extracted from plots like Fig.l) versus p for GOY and 
NS runs and their comparison with the SL formula Q and 
p/3. We obtain from our measured ^'p and the formula 
("p = ll("p/2 — 3p/2; this amplifies the error bars relative to 
Fig.l [inset(b)]. Error bars for ap are shown but not apparent 
since they are comparable to the symbol sizes. 



log{Hpq) = Dpq log(Gp) and \og{Hqp) = Dqp log(Gg); the 



scale factors 13, 
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Dqp are nonuniversal, but plots of 
log{Hpq) versus \og{Hqp) show data from all GOY and 
3d NS runs collapsing onto one universal curve within our 
error bars (Fig.0c for p = 6 and q = 9) for all k and Rex. 
Both ESS (Fig§) and GESS (Fig.l) remove the exponen- 
tial controlling factor [p^ from the leading asymptotic 
behavior of Sp in the far dissipation range and expose 
the remaining power-law dependence on k. Also, it is 
easy to see analytically that GESS plots (Fig.||) amplify 
slope differences between inertial- and dissipation-range 
asymptotes relative to ESS plots (Fig.^j). 

How robust is the fair agreement of (Fig-§) with 
the SL formula? Some studies p^-M suggest that, as 



Re\ — > 00, 5C,p = (p/3 — Qp) 0. Numerical solutions of 
the 3d NS equation can at best achieve 11 2^ Re\ ^ 
220, too small, by far, to resolve this issue, so we address 
it for the GOY model, by studying the range 4 x 10"* < 
Rex < 3 X 10^. We find (Fig.|) that 5C,p does not vanish 
with increasing Rey, if anything, it rises marginally [ pl| . 
Systematic experimental studies at high Rex are perhaps 
the best way to check if the trends of Fig.^ obtain in the 
3d NS case. 

We remark that, if we assume the hierarchy 
[Gp+i/Gp] = [Gp/Gp_i]T X [limp_ooGp+i/Gp]i-T with 
7^^ = 2/3 (whose real-space analog is equivalent ||l[ to 
the SL moment hierarchy for the energy dissipation Q|) 
and use |2^ Gp{k) « Cpk^p, we get a difference equa- 
tion for /3p identical to the SL one (our /3p is their 



— rp/3). This, when solved with the boundary condi- 
tions /3o ^ = and limp^oo{f3p+i - Pp) = 2/9, 
yields the SL formula (via C,p = —f3p + p^^/S). How- 
ever, our GESS yields [Gp+i/Gp] « C^[Gp/Gp-i]'^'' with 
= iCp+i - Cp - l/3)/(Cp - Cp-i - 1/3)- Superficially, 
this might seem to violate the hierarchy assumed above, 
but it turns out to be consistent with our GESS form, 
if Tp = 7 - 2(1 - 7)/[9(Cp - Cp-i - Cs)], which is pre- 
cisely the SL difference equation. Of course, our GESS 
form can hold with Cp 7^ Cp ^ ; Fig- @ shows the quality of 
agreement between our measured Cp and Cp^- 

We use a pseudospectral method [Q for our numerical 
solution of the incompressible 3d NS equation. We force 
the first two fc-shells, use a box with side = tt and 
64'' modes. Our dissipation term — (j/ + VHk^)k^ allows 
for both viscosity v and hyperviscosity vh- For time in- 
tegration we use an Adams-Bashforth scheme (step size 
5t) [Q. Parameters for our 3d NS runs NSl-4 are given 
in Table |, where = Lb jv^ms is the box-size eddy- 
turnover time and Tav the averaging time, after initial 
transients have decayed over a period tj. We use Rex = 
Vrms\lv, where A = {^^ E{k)dkl ]^ k'^E(k)dkYl'^ , 
vrms = [(2/31%) E{k)dk]^/^ and E{k) - S2{k)k^. 
All Sp{k) are averaged over shells of radius k. Care must 
be exercised in choosing St and the forcing amplitude, 
otherwise there is a slow, but systematic, stretching of 
the data points along the asymptotes in Figs. |l| and ^ 
with increasing Tav (over the time scales of our \ow-Rex 
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TABLE I. Parameters v (viscosity), vh (hyperviscosity), Re\ (Taylor-microscale Reynolds number), (box-size 
eddy-turnover time), Tav (averaging time), rt (transient time) and kd (dissipation-scale wavenumber) for our 3d NS runs 
NSl-4 (A:„a:r = 64) and GOY-model runs Gl-8 (fc,„aa; = l^'^ki^)- The step size(5t) used is 0.02 for NSl-4, 10"" for Gl-4, and 
2 • 10"^ for G5-8. 



Run 
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VH 


Rex 




Tt/Te 




kmax 1 kd 


NSl 


5 ■ 


IQ-* 





~ 3.5 


~ 3 ■ 10* 


~ 1 


2 


~ 4 


NS2 


2 ■ 


IQ-* 





~ 8 


~ 3 ■ 10* 


~ 1 


~ 2.5 


~ 2.3 


NS3 


5 ■ 


10-* 


5 • lO-** 


~ 3.5 


-3-10* 


~ 1 


~ 1 


~ 6.5 


NS4 


5 ■ 


10-* 


10-^ 


~ 22 


~ 3 ■ 10^ 


~ 10 


~ 7 


~ 2 


Gl-4 


5- 10" 







4 ■ 10* - 3 ■ 10^ 


~ (1.5 - 2.0)10* 


~ 500 


~ 2500 


~ 2^ - 2^ 


G5-8 


5- 10' 


8 _ 10-9 





3.5 • lO'"' -3-10'' 


~ (0.7- 1)10^ 


~ 500 


~ 2500 


~ 2^ - 1 
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FIG. 3. Log-log (base 10) plots of Ge versus (a) G15 and (b) Gg illustrating our fc-space GESS; (c) 
the universal inertial- to dissipation-range crossover (see text). The line shows the SL, inertial-range 




H(i,g versus 7^9,6 showing 
prediction. 



NS runs). Fortunately, this hardly affects our exponents: 
any attendant systematic errors in Fig. ^ are certainly 
less than the random errors indicated. Also, the agree- 
ment between our GOY and NS runs confirms our results. 
Our GOY-model data are, of course, of much better qual- 
ity. Here Fourier components of the velocity are labeled 
by a discrete set of wave vectors fc„ = koq"-. The dy- 
namical variables are the complex, scalar velocities w„ 
for each shell n; w„ is affected directly only by the ve- 
locities in nearest and next-nearest shells. In spite of its 
simplicity, this model yields scaling properties |9-12| akin 
to experimental ones. The GOY-model equations are: 



— Vn — iCn 

at 



Vklvn 



fni 



(4) 



where v is the kinematic viscosity, /„ the external force 

on sheU n, C„ = {ak„Vn+lVn+2 + + 

ckn-2Vn~iVn-2)* , and o, 6, and c can be fixed upto a con- 
stant by demanding ||ll[, for i^, fn = 0, that: Vn ^ fcn 
be a stationary solution of Eq.(4); and the GOY-model 
kinetic energy and helicity be conserved. We adopt the 
conventional parameters 11| fco — 2^*, q — 2, a — 
1, b = c = —1/2 and use /„ = 5 • 10-'^(1 -I- i)(5„,i, i.e., 
we force the first shell [ESl. The GOY-model structure 



functions are 5*^ 



ues of Cp obtain [gTj if we use S 



-Cp 



HI; reliable val- 



(|3[u„U„+lV„+2 + 

iw„w„+i/4] since this removes an underlying 



3— cycle. We have used E„^p to obtain Fig.|| |^^, but 
Sn,p in Figsj^^ for consistency with 3d NS. We use an 
Adams-Bashforth scheme |l^ (step size 6t) to integrate 
Eq. (4). The average of the time scale associated with 
the smallest wavenumber, (jwilfci)-^, gives the "box-size" 
eddy turnover time. Table ^ lists other parameters for 
our 8 GOY-model runs Gl-8, for which we use (cf., [M) 

E{k) = 5„,2/fc„, A = 2VfcoEn^n.2/Enfc«^n,2r', 

and v„ns = [fco I]« '5'„,2/7r]^/^. This yields Rex 
as expected at large Rex- 

Experimental evidence for the slope change in the 
dissipation range in real-space analogs of Fig.|l| was 
given by Stolovitzky and Sreenivasan who postu- 

lated iSp ~ iSg in the dissipation range and suggested 
a'p ~ (C3p/2+p/2)/(C9/2 + 3/2). We have not been able to 
obtain a simple relation between our ap and their a'^ (un- 
like |l^ that between (p and (p) since Sp does not have 
a power-law dependence on k in the dissipation range. 
It would be very interesting to extend such experimen- 
tal studies to test the universality of dissipation-range 
asymptotics (e.g., in different flows) and the crossover 
suggested here. The universal multiscaling in the dis- 
sipation range that we have elucidated is a manifesta- 
tion of strongly intermittent (multifractal) dissipation 
which is believed to occur [|5| even at low Rex- We be- 
lieve that this multiscaling should extend far enough into 
the dissipation range before corrections set in because of 
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FIG. 4. Log-log plot(base 10) of 5("p versus the Tay- 
lor-microscale Reynolds number Re\ for our GOY runs (Gl-8) 
with p = 6, 8 . . . , 20 (from bottom to top). The dotted (p — 6) 
and dashed {p = 8) lines show the SL results [^. Error bars 
are shown but are often smaller than the symbol sizes. 



the breakdown of (a) the incompressibihty assumption 
(at large Mach numbers) and/or (b) hydrodynamics (at 
molecular length scales). Preliminary studies (26| yield 
similar phenomena in MHD turbulence. 
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